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Context. The contribution of the merging process to the early phase of galaxy assembly at z > 1, and in particular for the build-up of 
the red sequence, remains to be accurately assessed. 
Aims. We aim to measure the major merger rate of star-forming galaxies at 0.9 < z < 1.8, using close pairs identified from integral 
field spectroscopy (IFS). 

Methods. We use the velocity field maps obtained with SINFONI/VLT on the MASSIV sample, selected from the star-forming 
population in the VVDS. We identify physical pairs of galaxies from the measurement of the relative velocity and the projected 
separation (r p ) of the galaxies in the pair. Using the well constrained selection function of the MASSIV sample we derive the gas-rich 
j.^ ■ major merger fraction (luminosity ratio fi = L 2 /Li > 1/4), and, using merger time scales from cosmological simulations, the gas-rich 

major merger rate at a mean redshift up to z = 1.54. 

Results. We find a high gas-rich major merger fraction of 20.8+^ 5 8 2 %, 20.1^°% and 22.0^ 3 3 7 % for close pairs with r p < 20/r' 
kpc in redshift ranges z = [0.94, 1.06], [1.2, 1.5) and [1.5, 1.8), respectively. This translates into a gas-rich major merger rate of 
0. 1 16+j; 038 °y r ' ' °- 147 -o 037 °y r ' and °- 127 -o 042 °y r 1 at z = 1 .03, 1 .32 and 1 .54, respectively. Combining our results with previous 
studies at z < 1, the gas-rich major merger rate evolves as (1 + z)", with n = 3.95 ± 0.12, up to z = 1.5. From these results we infer 
that ~ 35% of the star-forming galaxies with stellar masses = 10 10 - 10 10 5 M e have undergone a major merger since z ~ 1.5. We 
develop a simple model which shows that, assuming that all gas-rich major mergers lead to early-type galaxies, the combined effect of 
gas-rich and dry mergers is able to explain most of the evolution in the number density of massive early-type galaxies since z ~ 1.5, 
with our measured gas-rich merger rate accounting for about two-thirds of this evolution. 

Conclusions. Merging of star-forming galaxies is frequent at around the peak in star formation activity. Our results show that gas-rich 
mergers provide an important contribution to the growth of massive galaxies since z ~ 1.5, particularly on the build-up of the red 
■ sequence. 

Key words, morphology — galaxies: evolution — galaxies: mergers 
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1. Introduction The stellar mass density in the red sequence has in- 
creased by a factor of ~ 10 in the 2.5 Gyr between z = 

^ ■ Understanding the mechanisms involved in the mass assembly of 2 and z = 1, but only by a factor of ~ 2 in the 

galaxies and their relative role over cosmic time is an important last 7 - 8 Gyr of c osmic history (e.g ., lArnouts et al] 120071: 

open topic in modern astrophysics. In particular, the evolution of |y er gani et alj 120081: lllbert et alj 12010) . Major mergers, the 

the red sequence, which includes passive galaxies dominated by mer ger of two galaxies with similar stellar masses, is an effi- 

old stellar populations and an early-type (E/S0) morphology, im- dent mechanism t o create new passive early-typ e galaxies (e.g. 

poses fundamental constraints to formation and evolution mod- | Naab et all [20061: IRothberg & Joseph! [2006allbt i Hopkins et all 

els - I2OO8!; IRothberg & Fischerl 12010b iBournaud et all 1201 lb . Thus. 

the knowledge of the merger rate at z > 1 is an important in- 

* This work is based mainly on observations collected at the put to estimate the relative contribution o f merging and cold-gas 

European Southern Observatory (ESO) Very Large Telescope, Paranal, accretion (e.g., iDekel & Birnboim 2006) in the early assembly 

Chile, as part of the Programs 179.A-0823, 177.A-0837, 78.A-0177, of galaxies, and in particular the role of merging in the build-up 

75.A-0318, and 70.A-9007. This work also benefits from data products f the red sequence, 
produced at TERAPIX and the Canadian Astronomy Data Centre as 

part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- The evolution of the merger rate since z ~ 1 is now well 

rative project of NRC and CNRS. constrained by direct observations. The early measurements us- 
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ing photometric pairs dPatton et alJll997t iLe Fevre et al.l 



2000) 



2003 



or post-merger morphological signatures ( ConseliceT 
Jogee et al. 2009), have been superseded by spectroscopic mea- 
surements confirming physical pairs from the redshift measure- 
ment of both components of a major merger with a luminos- 
ity/m a ss rat io yu > 1 /4 (e.g.. ILin et alj 120081: Ide Ravel et alj 
20091 1201 11), as well as for m inor mergers down to /v = 
1/10 dLopez-Saniuan et al]|201 ll) . With a parametrization of the 
merger rate evolution following oc (1 + zf, it is observed that 
the major merger rate evolution depends on the luminosity an d 
on the mass of the galaxy sample (e.g., Ide Ravel et al. 2009), 
massive galaxies with M* > 10 11 M G having a higher merger 
rate, but with little redshift evolution (n ~ — 2), while lower 
mass galaxies with M* = 10 9 - 10 11 M have a lower merging 
rate but with stronger redshift evolution (n ~ 3 - 4). This mass 
dependency seems to explain some of the apparent discrepancy 
of merger rate measurements made from observations targeting 
different mass samples. 

Beyond z ~ 1, direct measurements of the merger rate are 
still limited. Previous attempts to measure the major merger 
rate at z > 1 have focused on the id entification of merger 
remnants from morphological signatures (IConselice et al.l l2008. 

on the study of projected close 



1201 lbt iBlucket 



orphologic 

iDj2012), 



pairs (Rvan et al. 2008; Bluck et al. 2009; Williams et al. 201 1; 
iMan et al.ll2012tlMarmol-Oueralt6 et alJl2012|:lLaw et al.ll2012h. 
or on indirect estimations dCameron & Pettitt 2012t |Puech et alj 
120121) . These studies find a high merger rate to z ~ 2-3 but 
with a large scatter between different measurements. However, 
these results are up to now solely based on photometric mea- 
surements which are increasingly hard to correct for contamina- 
tion along the line of sight as redshift increases. Another com- 
plication stems from the morphological evolution of galaxies, 
with show more irregular morphologies at high redshifts, and 
a wavelength dependency with more multi-c omponent object s 
present when observed in the rest-frame UV dLaw et al.ll2007l) . 
with some of these components possibly related to strong star- 
forming regions rather than different dynamical components. 

To improve on this situation, it is necessary to obtain spec- 
troscopic confirmation of the physical nature of the photometric 
pairs at z > 1 . In the last years, NIR Integral Field Spectrographs 
(IFS), like SINFONI on the VLT or OSIRIS on the Keck, have 
opened the possibility for a systematic study of the dynami- 
cal field around high redshift galaxies in the optical rest-frame. 
Some examples are the MASSIVn (Mas s Assembly Survey with 
SINFONI in VVDS. IContini et all2012l) survey at 0.9 < z < 1.8, 
the SINS0 (Spe ctroscopic Imaging survey in the Near-infrared 
with SINFONI. iForster Schreiber et al.ll2009l) survey at z ~ 2, or 
the Keck-OSIRIS dLaw et al.ll2009h . the A MAZE (Assessing the 
Mass- Abundance redshift -Z- Evolution. iMaiolino et alj|2008l) 
and the L SD (Ly man-break ga laxies Stellar populations and 
Dynamics, Mannucci et al. 2009) surveys at z ~ 3. 

The MASSIV survey has been designed to target the peak 
of the star-formation rate at 0.9 < z < 1.8, filling the gap be- 
tween higher redshift (z ~ 2) IFS surveys with those at z < 1, 
e.g., IMAGES (In termediate MAss Galaxy Evolution Sequence, 
lYanget al.ll2008l) . The MASSIV survey has targeted 84 star- 
forming galaxies at 0.9 < z < 1.8 with SIN FONI , drawn 
from the VVDS3 (VIMOS VLT Deep Survey, ILe Fevre et al] 
120051) survey. MASSIV has been used as a unique opportunity 
to study in details the dynamical state of 0.9 < z < 1.8 galax- 



1 http://www.ast.obs-mip.fr/users/contini/MASSIV/ 

2 http://www.mpe.mpg.de/ forster/SINS/sins_nmfs.html 
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ies (Epinat et al. 2012), their metallicity gradients (Oueyrel et al. 

or the evolution of the fundamental mass-size-velocity re- 
lations since z ~ 1.2 dVergani et al.ll2012l) . 

In this paper, using the MASSIV survey, we present for the 
fist time a measurement of the gas-rich major merger rate of star- 
forming galaxies from kinematical close pairs at 0.9 < z < 1.8. 
Thanks to the large field-of-view of IFS we have access to the 
complete surrounding volume of the galaxies when searching 
for close kinematical companions. In addition, the well-defined 
selection of sources from the VVDS and the well controlled 
selection function of MASSIV observations ensures the study 
of a represent ative population of star-forming galaxies at these 
redshifts (see IContini et al.l|2012l for details). This all together 
enables the measurement of average volume quantities like the 
merger fraction and rate. 

The paper is organized as follows: in Sect. [2] we summa- 
rize the MASSIV data set used to identify merging pairs, and in 
Sect. [3] we develop the methodology to measure the merger frac- 
tion from IFS data. We report the gas-rich major merger fraction 
in MASSIV in Sect. [4] and derive the gas-rich major merger rate 
in Sect. [5] We discuss the implication of our results in Sect. [6] 
Finally, we present our conclusions in Sect. [7] We use Hq = 100k 
km s _1 Mpc~', h = 0.7, Q m = 0.3, and Q A = 0.7 throughout this 
paper. Al l magnitudes ref er to the AB system. The stellar masses 
assume a SalDeterl (119551) initial mass function (IMF). 

2. The MASSIV data set 

The galaxy sample studied in this paper is the final release of 
the MASSIV project (ESO Large Programme 179.A-0823; PL: 
T Contini). A full description of the sample can be found in 
IContini et al.l (120121) . We briefly summarise some properties of 
this sample of 84 galaxies below. The galaxies were selected 
from the VVDS in t he RA = 2h area of the deep (/ A b < 24; 
ILe Fevre et al.l 120051) and ultradeep (7ab ^ 24.75; Le Fevre et 
al., in prep.) surveys, and in t he RA = 14h and R A = 22h of 
the wide survey (7 A b < 22.5: iGarilli et al.l 12008b . The galax- 
ies were chosen according to their redshift such that their Ha 
line (or [O Iii]/15007 in a few cases) was visible in the J or H 
band, and was not affected by a bright OH sky-line. Galaxies 
were randomly selected to be star-forming on the basis of their 
[O ii]/13727 emission line strength and equivalent width at z < 
1.5 (Fig. [H>, and of their re st-frame UV spectral properties at 
z > 1.5 dContini et al.ll2012l) . The observations have been per- 
formed between April 2007 and January 2011. Most (85%) of 
the galaxies in the sample have been observed in a seeing-limited 
mode (with a spatial sampling of 0.125 arcsec/pixel). However, 
eleven galaxies have been acquired with Adaptive Optics as- 
sisted with a Laser Guide Star (AO/LGS, seven with 0.05 and 
four with 0.125 arcsec/pixel spatial sampling). The data reduc- 
tion was performed with the ESO SINFONI pipeline (version 
2.0.0), using the standard master calibration files provided by 
ESO. The absolute astrometry for the SINFONI data cubes was 
derived from nearby bright stars also used for point spread func- 
tion (PSF) measurements. Custom IDL and Python scripts have 
been used to flux calibrate, align, and combine all the individual 
exposures. For each galaxy a non sky-subtracted cube was also 
created, mainly to estimate the effective s pectral resolu t ion. F or 
more details on data reduction, we refer to Epinat et aTJ d2012l) 

We use the Ha emission line (or [O iii]/15007 in a few cases) 
to derive the kinematical maps (flux, velocity field and veloc- 
ity dispersion map) of the MASSIV galaxies. To estimate their 
dynamical properties, we assume that the ionized gas rotates in 
a thin disc with two regimes for the rotation velocity, a solid 
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Fig. 1. MASSIV selection in the [On]/l3727 flux vs equivalent width plane (see IContini et al.ll2012l for details). Red squares are 
the VVDS sources at 0.94 < z < 1.5 with individual [Oll]/l3727 measurements in the VVDS-Wide 14h (top-left), VVDS-Wide 
22h (top-right), VVDS-Deep (bottom-left), and VVDS-Ultradeep surveys (bottom-right). The dashed lines mark the selection of 
MASSIV star-forming galaxies. White dots are those VVDS sources that fulfil the MASSIV selection. Green pentagons are the 
MASSIV galaxies observed with SINFONI/VLT. [A colour version of this plot is available at the electronic edition]. 



body shape in the innermost regions and a plateau in the out- 
skirts. Using ax 2 minimization we produce seeing-corrected ve- 
locity and dispersion maps of galaxies with geometrical inputs 
weighted for the signal-to-noise ratio of each pixel. We esti- 
mate the geometrical parameters used in the fitting model on the 
t-band best-seeing CF HTLS Megacam images for all galaxies 
dGoranova et al J2009l) . except for VVDS-Wide 1 4h galaxies that 
were covered with the CFHT-12K/CF HT camera jLe Fevre et all 
120041) . We use the GALFIT software dPeng et alj|2002l) that coi> 
volves a PSF with a model galaxy image based on the initial pa- 
rameter estimates fitting a Sersic ( 1968) profile. Residual maps 
from the fitting were used to optimize the results. At the end 
of the fitting procedure, GALFIT converges into a final set of 
parameters such as the centre, the position angle, and the axial 
ratio. The /-band images were also used to correct for SINFONI 
astrometry, using the relative position of the PSF star. The mor- 
phology and kinematics maps o f the 50 first-epoch galaxies are 
presented in lEpinat et all d2012l) together with a more extensive 
discussion on the model fitting procedure. The 34 second-epoch 
galaxies, already included in the present work, will be presented 
in a future paper. 

We obtain the stellar mass of MASSIV galaxies from a spec- 
tral energy distribut ion (SEP) fit to the phot ometric and spectro- 
scopic data with iBruzual & Charlotl d2003l BC03) stellar pop- 
ulation synthesis models using the GO SSIP2 Spectra l Energy 
Distrib ution tool dFranzetti et al.1 12008). We assume a[Salpeter 
d 19551) IMF, and a set of delayed exponential star formation his- 



tories with galaxy ages in the range from 0. 1 to 15 Gyr. As input 
for the SED fitting, in addition to the VVDS spectra, we use 
the multi-b and photometric ob servations available in the VVDS 
fields (see IContini et all [20121 for further details). Following 
Walc heretalT ?2008) we adopt the probability distribution func- 
tion to obtain the stellar mass. 

In all, the final 84 MASSIV galaxies are representative of the 
star-forming (star formation rate S FR > 5M yr _1 ) population 
of = 10 9 - 10 11 M Q galaxies at 0.9 <z< 1.8. 



3. Measuring the merger fraction from IFS data 

We define as a close pair two galaxies with a projected sepa- 
ration in the sky plane r™ n < r p < r™ ax and a rest-frame rel- 
ative velocity along the line of sight Av < Av max . We used 
Av max = 500 kms- 1 and r™ x = 20-30/r 1 kpc (see Sect.E2 for 

details), while setting r™ n = 0. We therefore searched for close 
companions in the kinematical maps of the MASSIV sources, 
analysing and classifying the sam ple using the veloci ty field and 
the velocity dispersion map (see Epinatet a~D l2012l for details 
about the classification). We find 20 close pair candidates in the 
MASSIV data cubes, and we study these systems in details to 
split them into major (luminosity difference between both com- 
ponents n = La/L] > 1 1 A) and minor (/i < 1/4) close pairs 
(Sect®. 
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If there was N p major close pairs in our sample, the gas-rich 
major merger fraction is 



xlO" 



f - Np 
/mm - — , 



(1) 



where N is the number of principal galaxies targeted for the sur- 
vey. We named principal galaxy the source in the pair closest to 
the kinematical centre of the targeted system, even if it is not the 
brightest/more massive galaxy in the pair. The simple definition 
in Eq. (Q]l is only valid for volume-limited samples. While our 
sample is not only luminosity-limited but spectroscopically de- 
fined, we must take into account the different selection effects, 
both in the VVDS parent samples and in MASSIV, in our com- 
putation of the merger fraction. 

3.1. Accounting for selection effects in the VVDS 

Since a fraction of the total number of potential targets in the 
VVDS fields have been spectroscopically observed and since 
the redshifts are not measured with 100% accuracy, we must 
correct for the VVDS Target Sampling Rate (TSR) and the 
Redshift Success Rate (SSR), computed as a function of redshift 
and source magnitude. The SSR has been assum ed indepen- 
dent o f galaxy type, as demonstrated up to z ~ 1 in lZucca et al.l 
(2006). E very VVD S source has a redshift confidence flag (see 
iLe Fevre et al.ll2005l for details), that can be flag = 4 (redshift 
99% secure), flag = 3 (95% secure), flag = 2 (75% secure), flag 
= 9 (redshift from a single emission line), flag = 1 (50% se- 
cure), or flag = (no redshift information). As several VVDS- 
Deep galaxies with flag = 2 have been re-observed in the VVDS- 
Ultradeep survey, providing a robust measurement of their red- 
shift, this offers the opportunity to estimate the reliability of 
VVDS-Deep flag = 2 sources. We thus define a weight W29 to 
take this into account. We also define the weight W29 for flag = 
9 sources by comparison with the latest photometr ic redshifts 
in the VVDS-Deep field (see ICucciati et ail 1201 PL for details 
about the latest photometric data set in the 2h field). By defi- 
nition, the W29 weight is equal to 1 for flag = 3 and 4 sources 
in VVDS-Deep, and for all sources in VVDS-Wide and VVDS- 
Ultradeep. We derived the spectroscopic completeness weight 
for each galaxy i in the VVDS catalogue as 



"VVDS 



(z, Tab, flag) 



w 29 (z,/ A B,flag) 



T SR'(z, 1 ab) xSSRKz, 1 ab) 



(2) 



The TSR, SSR and w 29 on VVDS-D eep and VVDS- 
Ultradeep were mea sured in previous works (lllbert et al.ll2006l 
ICucciati et ail 120 121). We assume TSR = 0.22 in VVDS-Wide 
fields dGarilli e t al. 2008), and we detail the computation of the 
S S R in VVDS-Wide fields in Appendix|A] 

3.2. Accounting for selection effects in the MASSIV survey 

The MASSIV sources were drawn from the star-forming galaxy 
population in the VVDS. We correct for three basic selection 
effects in MASSIV. 

- The selection weight, w se i. We define this weight as the frac- 
tion of star-forming galaxies that fulfil the MASSIV selec- 
tion in the log/[ 0n ] [10~ 17 erg s -1 cirT 2 ] vs log(-£W[0 11]) 
[A] plane at z < 1 .5 (Fig.|TJ. The weight w se i tells us how rep- 
resentative the MASSIV selection of the global star-forming 
population in each of the VVDS surveys (Wide, Deep and 
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Fig. 2. Typical field-of-view of a MASSIV source. We show the 
flux of the source 220397579 in the five channels of its reduced, 
flux calibrated, data cube centred on the position of the Ha emis- 

sion line (scale in the right). The pixel spatial scale is 0.125 , 
that at the redshift of the source (z = 1 .0379) corresponds to ~ 1 
kpc. The source is at the centre of the image (black cross), while 
the white circle delimits the 20/z~' kpc area around the source. 
In this particular case, 99% of this area is covered with MASSIV 
data and the two negatives of the source (black regions) are also 
excluded. 



Ultradeep) is, and gives more importance to the more rep- 
resentative samples. The VVDS-Deep and Ultradeep have 
w se i ~ 0.67 at z < 1.5, while the VVDS-Wide fields have 
w S ei ~ 0.53 in the same redshift range. We assume w se i = 1 
at z > 1.5, where the selection is based on colour/spectral 
properties and all the star-forming galaxies in the VVDS are 
thus pre-selected; 

the MASSIV IFS Rate (MIR) is defined as the fraction of 
galaxies that fulfil the MASSIV selection and which were 
finally observed with SINFONI (Fig. Q}. The MIR ranges 
from 0.43 for VVDS-Ultradeep to 0.05 for VVDS-Deep; 
the MASSIV Success Rate (MSR) is the fraction of observed 
sources with a reliable kinematical classification. This frac- 
tion is always high, MSR > 0.8. 



Finally, the MASSIV weight is 

Wsel (X, Z) 



"MASSIV 



(x,z) 



MIRJ(x,z)xMSRj(x,z) 



(3) 



where the index j spans for the MASSIV sources and x refers to 
the VVDS survey (Wide at 14h, Wide at 22h, Deep or Ultradeep) 
to which the source belongs. 

3.3. Area correction 

Thanks to the large field-of-view of SINFONI, we have access to 
the complete surrounding volume of the principal galaxy when 
searching for close companions. This is a major advantage with 
respect to long-slit spectroscopic surveys, in which the observed 
number of close pairs is proportional to TSR 2 , diminishing the 
statistics. However, the final reduced MASSIV data cubes cover 
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Table 1. Close pair candidates at 0.9 < z < 1.8 in the MASSIV sample. 



ID 


z 


'"n 


Av 


Am, 


Ara^ 


Classification 






(n Kpc) 


(Km s ) 








020294045 


1.0028 


2.9 


180 


0.4 


0.3 


Major merger 


020386743 


1.0487 


3.7 


80 


1.8 




Minor merger 


020461235 


1.0349 


2.8 


15 


1.7 




Minor merger 


140096645 


0.9655 






> 1.5 




Minor merger 


220397579 


1.0379 


14.4 


340 


0.4 


-1.4 


Major merger 


220544394 


1.0101 


7.1 


50 


1.3 




Major merger 


020167131 


1.2246 


15.2 


130 


0.2 


0.1 


Major merger 


020218856 


1.3103 






> 1.5 




Minor merger 


020240675 


1.3270 






> 1.5 




Minor merger 


020283083 


1.2818 


3.8 


5 


0.7 




Major merger 


020283830 


1.3949 


8.5 


500 


1.9 




Minor merger 


020465775 


1.3583 


3.6 


40 


0.7 




Major merger 


220376206 


1.2445 


13.4 


400 


2.4 




Minor merger 


220544103 


1.3973 


5.0 


75 


-1.1 




Major merger 


910154631 


1.3347 


4.2 


130 


0.8 




Major merger 


910296626 


1.3558 


12.1 


165 


-0.1 


-0.2 


Major merger 


910337228 


1.3955 


9.5 


220 


1.4 




Major merger 


020116027 


1.5302 


26.8 


100 


0.7 


0.5 


Major merger 


910186191 


1.5399 


12.7 


450 


-0.2 


-2.4 


Major merger 


910274060 


1.5694 


3.4 


10 


0.2 




Major merger 



a finite area in the sky plane. To deal with this, we define the area 
weight as 



A r 



^MASSIV 



(4) 



where A r is the area subtended in the sky plane by a circle of 
radius r™ ax at the redshift of the source j, and Amassiv is the 
area of the same circle covered by the reduced SINFONI mosaic 
(Fig. [2]). We state that the optimum search radius is r™ x = 20h~ l 
kpc. This choice of r™ x minimizes the area correction, with 

w area ~ 1 m a U cases. In addition, the two negatives of the main 
source produced by the offset observing procedure and that ap- 
pear in the extremes of the reduced SINFONI mosaic are then 
excluded from the search area. 

Finally, the corrected gas-rich major merger fraction is 



7mm 



k ? W VVDS W MASSIV W area 



Z5X 



(5) 



VVDS ^MASSIV 



where the index j and k spans respectively for all MASSIV 
galaxies and for the MASSIV galaxies with a major close com- 
panion. The error budget in the major merger fraction is domi- 
nated by t he low statistics. We use the Bayesian approach from 
Cam eronl d201 lb to measure the statistical error in the raw ma- 
jor merger fraction, i.e., in N p /N, then scale it with the previous 
weighting scheme. 

4. Mergers classification and the gas-rich major 
merger fraction in MASSIV 

In this section we measure, for the first time, the major merger 
fraction at 0.9 < z < 1.8 from spectroscopically-confirmed close 
pairs. The MASSIV observational strategy defines three natural 
redshift bins. The low redshift MASSIV sources were observed 
in the J band, while the higher redshift ones in the H band. 
This translates to a gap in the redshift distribution at z ~ 1.1. 



In addition, the selection function of MASSIV targets changes 
at z = 1.5 (Sect. 0, providing another redshift boundary. We 
take advantage of these natural separations in the data to esti- 
mate the gas-rich major merger fraction in three redshift ranges, 
z r ,i = [0.94, 1.06], z r ,2 = [1.2,1.5), and z r3 = [1.5,1.8). We 
restrict our study to those galaxies with /ab ^ 23.9 to ensure 
completeness in the detection of close pairs (see Appendix IB~1 
for details). 

We follow the steps described bellow to split the close pairs 
candidates in the MASSIV data cubes into major and minor 
mergers: 



1. 



As described in Ep inat et al.l d2012l) . we have performed a 
classification of the MASSIV sources based on the shape of 
the velocity field (regular or irregular) and the close environ- 
ment (isolated or not-isolated). From this classification, we 
have pre-selected as close pair candidates the non-isolated 
sources a nd those identifie d as mergers from the velocity 
field (see lEpinatet~aT]l2012l for more details). We identified 
20 close pair candidates (Table[T]i. 

To study in details these close pairs candidates we used the 
deepest /-band images in the VVDS fields: CFTH12K (14h - 
field, exposure time of f exp = 3.6 ks. iLe Fevre et al. Il2004l) 
CFHT LS-Wide (22h field, f e x P = 4 - 10 ks. lGoranova et all 
1 20091). and CFTHLS -Deep (02h field, f exp ~ 300 ks, 
iGoranova et alJ l2009l) . We run SExtractor on the systems 
with well separated sources, and GALFIT (with two Seisic 
components) on the blended ones, to estimate the luminos- 
ity difference in the i band between both sources, Am, = 
mi,2 - We take Am, < 1.5 (factor four or less in lu- 
minosity) to identify major mergers and Am,- > 1 .5 for minor 
mergers. The observed i band corresponds to ~ 300-350 nm 
rest-frame in the redshift range of our sample. W e stress that 
for th e blended sources, we run GALFIT v3.0 dPeng et al.l 
2010) without imposing any constraint to the parameters of 
the fit and we only use the information from the kinematical 
maps to set the initial positions of the sources. 



5 



C. Lopez-Sanjuan, et al.: MASSIV: the major merger rate at 0.9 < z < 1.8 from IFS-based close pairs 



VVDS-Deep: 020294045 
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Fig. 3. The /-band image {left panel, scale in ADU in the right), the velocity field {central panel, scale in km s in the right), 
and the velocity dispersion map {right panel, scale in km s in the right) of the MASSIV source 020294045 (major merger). On 
each map, North is up and East is left. The level contours mark the isophotes of the two components obtained with GALFIT in the 
/-band image. The principal galaxy (red/white) is the one closer to the kinematical centre of the system and set the origin in right 
ascension {a) and declination (<5), while the companion (blue/black) is the secondary component. The outer contour marks the 3.43 
ADU (7cr SKy ) isophote. The next contours mark brighter isophotes in 1.96 ADU (4cr sky ) steps. The luminosity difference between 
both components in the i band, Am,, as well as the CFHTLS exposure time, f exp , are shown in the left panel. The redshift of the 
source and the on-source SINFONI exposure time, /sinfoni, are shown in the central panel. [A colour version of this plot is available 
at the electronic edition]. 
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Fig. 4. The same as Fig. [3] but for the ^T s -band image {left panel) of the MASSIV source 020294045 (major merger). The outer 
contour marks the 6.76 ADU (2<x SKy ) isophote, while brighter isophotes increase in 5.07 ADU (1.5cr sky ) steps. The flux difference 
between both components in the K s band, Am^, as well as the WIRDS exposure time, f exp , are show in the left panel. [A colour 
version of this plot is available at the electronic edition]. 



3. We confronted the images and the two component fits from 
GALFIT with the velocity field and the velocity dispersion 
map of the sources. We compared the distribution of Ha 
emission and the geometry of the velocity field to the rest- 
frame UV continuum (or rest-frame visible when NIR im- 
ages are available). The presence of two components with 
position and geometry concordant in the velocity field and in 
the continuum images, is a strong indication of the reality of 
the pair. 

4. We run SExtractor in the residual image from GALFIT with 
the principal source subtracted to obtain a second estima- 
tion of nij i, while with the companion source subtracted to 
estimate my. Then we compared the Am, derived from the 
GALFIT modeling with that from these SExtractor estima- 
tions. We find good agreement between both measurements 
(difference of ~0.2 mag or less). The major/minor merger 
classification does not change from the initial estimate. 

5. Finally, we also explored Am# s for well separated galax- 
ies and for one blended source (020294045). We used the 
K s -band image s from UKIDSS-DXS survey (22h field , 
iDveet al.i r2006) and WIRDS (2h field. iBielbv eta!ll20TTh . 



The observed K s band corresponds to ~ 0.8 - 1 yum rest- 
frame in the redshift range of our sample, and is a better 
tracer of the stellar mass content of the galaxies. After this 
second check there is not change in the major/minor merger 
classification, except for one close pair (source 910186191), 
supporting the previous /-band results. 

In the next sections we present the decomposition and the 
major/minor classification of the 20 close pair candidates in the 
MASSIV data cubes (Tableffl). 

4.1. Close pair candidates at 0.94 < z < 1.06 

The weighted mean redshift of the first redshift bin is z7j = 1.03. 
We have identified 6 close pairs over 1 8 sources in this redshift 
range: 

- 020294045 (Fig. [3). Major merger. The velocity map sug- 
gests two projected components. The companion is toward 
the north and presents a steep velocity gradient compared 
with the principal galaxy. The GALFIT model with two com- 
ponents in the i band reproduces both the position of the two 
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VVDS-Deep: 020386743 
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Fig. 5. The same as Fig. [3] but for the MASSIV source 020386743 (minor merger). The outer contour marks the 1 .41 ADU (2<x s ky) 
isophote, while brighter isophotes increase in 1.76 ADU (2.5<x sky ) steps up to 10 ADU to avoid crowded figures. [A colour version 
of this plot is available at the electronic edition]. 



VVDS-Deep: 020461235 
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Fig. 6. The same as Fig. [3] but for the MASSIV source 020461235 (minor merger). The outer contour marks the 1.31 ADU (3o" sky ) 
isophote, while brighter isophotes increase in 1.75 ADU (4cr sky ) steps up to 7 ADU to avoid crowded figures. [A colour version of 
this plot is available at the electronic edition]. 



VVDS-Wide: 140096645 
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Fig. 7. Left and central panels are the same as in Fig. [3] but for the MASSIV source 140096645 (minor merger). The level contours 
mark the isophotes of the original /-band image. The Ho- companion source is not detected in the i band. The outer contour marks the 
2.64 ADU (lcr^y) isophote, while brighter isophotes increase in 7.94 ADU (3cr sky ) steps up to 45 ADU to avoid crowded figures. 
The right panel shows the detection fraction of fake galaxies injected in the /-band image as a function of the magnitude of the 
companion source m,-.2- The size of the dots marks the detection curves for extended, normal and compact galaxies (see text for 
details). The vertical line marks the limiting magnitude for major companion (Am, = 1.5). [A colour version of this plot is available 
at the electronic edition]. 



kinematical components in the velocity map and the high 
velocity dispersion in the overlapping region between both 
components. The luminosity difference is Am, = 0.4. The 
separation between the sources is 2.9h~ l kpc and their rel- 
ative velocity is Av ~ 180 km s _1 . Despite the close sep- 
aration lead into significant overlap, this system also gives 
a satisfactory fit in the WIRDS A^-band image (Fig]5j>. In 



this case Aiuk s = 0.3, confirming that this system is a major 
merger. 

- 020386743 (Fig. O. Minor merger. The velocity map sug- 
gests two projected components. The companion is toward 
the north and presents a different velocity than the northern 
part of the principal galaxy. The GALFIT model with two 
components suggests that we are only detecting the western 



7 



C. Lopez-Sanjuan, et al.: MASSIV: the major merger rate at 0.9 < z < 1.8 from IFS-based close pairs 



< 




1.5 1.0 0.5 0.0 -0.5 -1.0 -1.6 -2.0 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 -2.0 1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 -2.0 

A« (") Ac (") A a (") 



Fig. 8. The same as Fig. [3] but for the MASSIV source 220397579 (major merger). The level contours mark the isophotes of the 
original /-band image. The outer contour marks the 5.36 ADU (3oSky) isophote, while brighter isophotes increase in 2.68 ADU 
(1.50",^) steps. [A colour version of this plot is available at the electronic edition]. 
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Fig. 9. The same as Fig. [3] but for the MASSIV source 220544394 (major merger). The outer contour marks the 4.3 1 ADU (2o~ s k y ) 
isophote, while brighter isophotes increase in 3.23 ADU (1.5cr S ky) steps. [A colour version of this plot is available at the electronic 
edition]. 



part of the companion galaxy, with the eastern part being 
too faint. This is also consistent with the velocity dispersion 
map, which shows a regular pattern in the overlapping re- 
gion. The separation between the sources is 3.7 hr l kpc and 
the luminosity difference is Am, = 1.8. Hence, we classify 
the system as a minor merger. 

- 020461235 (Fig.|6]i. Minor merger. The velocity map sug- 
gests two projected components. The companion is toward 
the south-west and presents a nearly constant velocity, in 
contrast with the velocity gradient of the principal galaxy. 
The GALFIT model with two components suggests that the 
companion is highly distorted, perhaps because the system is 
in an advanced merger stage (i.e., after first pericenter pas- 
sage or pre-coalescence). The principal galaxy also presents 
slightly distorted isophotes in the direction of the compan- 
ion. The luminosity difference from the GALFIT model is 
Am, = 1.7, while from the residual maps is Am,- = 1.9. This 
difference, in combination with the very different levels of 
distortion in the principal and the companion, strongly sug- 
gests that this system is a minor merger. The separation be- 
tween the components is 2.8/i _1 kpc. 

- 140096645 (Fig. |7). Minor merger. A small companion is 
identified in the velocity map about 1.2 arcsec away from 
the principal galaxy. There is no continuum detection in the 
/-band image in the position of the Ha companion, although 
there is a ~ lcr sky excess emission. To estimate if the non- 
detected companion could be bright/massive enough to lead 



to a major merger, we injected fake companion sources in the 
/-band image and studied their detection fraction as a func- 
tion of the fake companion luminosity m,^. The fake sources 
were modelled with a Sersic function and convolved with a 
typical PSF of the 14h field using GALFIT. For each fake 
source we assumed a random inclination and position angle, 
a Sersic index n s = 1, i.e., an exponential disc (we checked 
that the detection curve is similar assuming either n s = 0.5 
or n s = 2), and an effective radius r e given by the r e - M+ 
relation i n MASSIV, log r R = 0.36 + O.37[log(M*/M ) - 10] 
(see also lVergani et al.|2012l) . To estimate the stellar mass of 
the fake companion, we took Am,- as a proxy of the mass ra- 
tio between the principal galaxy, for which the stellar mass 
is known, and the companion. Then, we applied Poissonian 
noise to the model and injected it in the expected posi- 
tion of the possible companion. We measured, for different 
luminosities m, 2, which fraction of the 500 injected fake 
sources were detected. We repeated the previous steps for 
fake sources with 2r e (extended sources) and 0.5r e (compact 
sources), spanning all the possible sizes of the real sources. 
We show the result of this experiment in the right panel of 
Fig. 13 

We find that even for extended sources, we are ~ 90% com- 
plete at m,- 2 ~ 24, while the limiting magnitude for a major 
companion is mJ^ M = 23.8 (vertical line in the right panel 
of Fig. 17). Because of the small probability of non detection 
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Fig. 10. The same as Fig. [8] but for the MASSIV source 020167131 (major merger). The outer contour marks the 3.43 ADU (10<x sky 
isophote, while brighter isophotes increase in 2.74 ADU (8cr sky ) steps. [A colour version of this plot is available at the electronic 
edition]. 
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Fig. 11. The same as Fig. [7] but for the MASSIV source 020218856 (minor merger). The Ho- companion galaxy is not detected in 
the i band. The outer contour marks the 0.34 ADU (lcr sky ) isophote, while brighter isophotes increase in 1.01 ADU (3<x sky ) steps. 
[A colour version of this plot is available at the electronic edition]. 
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Fig. 12. The same as Fig. [7] but for the MASSIV source 020240675 (minor merger). The Ha companion galaxy is not detected in 
the i band. The outer contour marks the 0.36 ADU (l<x sky ) isophote, while brighter isophotes increase in 1.08 ADU (3<x sky ) steps up 
to 8 ADU to avoid crowded figures. [A colour version of this plot is available at the electronic edition]. 



of a major companion, we classify this system as a minor 
merger. 

- 220397579 (Fig. [8]). Major merger. The velocity map shows 
two different components with a separation of r p - 14.4/i~' 
kpc and a relative velocity of Av ~ 340 km s _1 . The com- 
panion is toward the north-west. The luminosity difference 
is Am, = 0.4, suggesting a major merger. The difference in 
the K s band is Am# s = -1.4. The negative sign implies that 
the companion is more luminous than the principal, which is 
brighter in the NUV rest-frame and in Ho-. This is consistent 
with the measured integrated metallicity of these sources, 



that is higher for the companion (Oue vrel et al.ll2Q12h . The 
suggested picture is that the system comprises a nearly face- 
on principal galaxy with intense star formation and low dust 
reddening, and a nearly edge-on companion galaxy with ei- 
ther a low level of star formation or strong dust reddening. 
We classify the system as a major merger. 
- 220544394 (Fig.|9]i. Major merger. The velocity map shows 
two different components with a separation of r p = 7. l/i 1 
kpc and a relative velocity of Av ~ 50 km s _1 , that are also 
well recovered by the GALFIT model with two components. 
The companion is toward the north, while the luminosity dif- 
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ference is Am,- = 1.3. We classify the system as a major 
merger. 

In summary for this redshift range, we classify 3 of the 6 
close pairs candidates as major mergers, and the three other ones 
as minor mergers. This translates to a gas-rich major merger 
fraction of /mm = 0.208+°;^ a t Tr,j = 1.03. 

4.2. Close pair candidates at 1.2 < z < 1.5 

The weighted mean redshift of the second redshift bin is z\- ,2 = 
1.32. This is a redshift range where there is no measurement 
of the major merger fraction from spectroscopic close pairs yet. 
In this framework, MASSIV provides an unique opportunity to 
measure the major merger fraction at this crucial epoch of galaxy 
evolution. We identify 1 1 close pair candidates over 30 galaxies 
in this redshift bin: 

- 020167131 (Fig. [TOl , Major merger. The velocity map 
shows two different components with a separation of r p = 
\5.2h~ l kpc and a relative velocity of Av ~ 130 km s . 
The line targeted in this case is [O iii]/15007, which explains 
the low signal. The companion galaxy is toward the south- 
east, while the luminosity differences are Am, = 0.2 and 
Am# s =0.1. We classify the system as a major merger. 

- 020218856 (Fig. [TT). Minor merger. A faint companion is 
detected in the velocity map about 1 .6 arcsec from the princi- 
pal galaxy. There is no continuum detection of this compan- 
ion at the depth of the /-band image. We followed the same 
steps than for the source 140096645 to estimate the detec- 
tion probability of companion galaxies, finding that we are 
100% complete up to /«, 2 = 26. We note that the curves for 
extended, normal and compact galaxies are similar, reflect- 
ing that the size of the PSF is larger than the assumed size 
of the fake sources. Because a major companion should be 
brighter than m^, M = 25.32 (vertical line in the left panel of 
Fig. rm>. and the non detection in the continuum means that 
the companion should have Am* > 2, we classify the system 
as a minor merger. 

- 020240675 (Fig. [T2l . Minor merger. A faint companion is 
detected in the velocity map about 2.4 arcsec from the princi- 
pal galaxy. There is no continuum detection of this compan- 
ion at the depth of the /-band image. We followed the same 
steps than for the source 020218856 to estimate the detection 
curves of companion galaxies. We find similar completeness 
curves. Because a major companion should be brighter than 
m^, M = 25.15 (vertical line in the left panel of Fig.[T2l. and 
the non detection in the continuum means that the compan- 
ion should have Am, > 2.5, we classify the system as a minor 
merger. 

- 020283083 (Fig.Qj). Major merger. The velocity map sug- 
gests two projected components separated by 3.8/z _I kpc and 
~ 5 km s _1 , with an extended region in the north-west. The 
GALFIT model with two components reproduces the shape 
of the velocity map and suggests that the extended region is 
due to the companion galaxy. The luminosity difference from 
the models is Am,- = 0.7. Thus, we classify this system as a 
major merger. 

- 020283830 (Fig. Q3). Minor merger. The velocity map 
shows two different components with a separation of r p = 
8.5/T 1 kpc and a relative velocity of Av ~ 500 km s . The 
companion only presents six detected pixels in the Ha map, 
already suggesting a minor companion. We used GALFIT to 
model the system, finding Am, = 1 .9. However, the model of 



the companion provided by GALFIT is a point-like source, 
thus overestimating the luminosity of the extended compan- 
ion. The measurement with SExtractor in the image with the 
principal galaxy subtracted suggests Am,- = 2.1. Since the 
companion is fainter than the major merger limit, we clas- 
sify the system as a minor merger. 

- 020465775 (Fig.[T3Tl. Major merger. The velocity map sug- 
gests two projected components. The GALFIT model with 
two components finds that the companion galaxy is toward 
the north-west of the principal galaxy, at r p = 3.6/i~' kpc and 
Av ~ 40 km s -1 . The Ho- emission of this system is located 
in the central part of the galaxies, since the /-band models 
match the velocity map at 90"^ level, i.e., we do not detect 
emission from the outer parts of the galax. The position of 
the companion explains the abnormal velocity pattern and 
the high velocity dispersion peak in the maps. The luminos- 
ity difference of the system is Am,- = 0.7, so we classify it as 
a major merger. 

- 220376206 (Fig. [TBI . Minor merger. The velocity map 
shows two different components with a separation of r p = 
13.4/! -1 kpc and a relative velocity of Av ~ 400 km s -1 . The 
companion, located toward the north, presents ten detected 
pixels and its velocity is inconsistent with that expected from 
the velocity field of the principal galaxy. We find Am,- = 2.4, 
so the system is a minor merger. 

- 220544103 (Fig.[J7]i. Major merger. The velocity map sug- 
gests two projected components. The southern component 
presents a large velocity gradient and defines the kinematical 
centre of the system, while the northern component is more 
extended and has a nearly flat velocity field. The GALFIT 
model with two components recovers the configuration in 
the velocity map and suggests that the southern component 
is edge-on, while the northern component is nearly face-on, 
thus explaining the observed high and null velocity gradi- 
ents. The separation between the components is r p = 5h~ l 
kpc, while their relative velocity is Av ~ 75 km s _1 . From 
the GALFIT models we estimate Am, = -1.1, and we clas- 
sify the system as a major merger. 

- 910154631 (Fig. [18]!. Major merger. The velocity map sug- 
gests two projected components, with the companion toward 
the north-west. The /-band image shows another two well 
separated sources close to the MASSIV target. These sources 
are not detected in Ha. To avoid contamination from these 
sources in the /-band photometry, we performed a four com- 
ponent fitting, with two components for the MASSIV target 
and one component for each nearby source. The GALFIT 
model finds the second component of the MASSIV target at 
the expected position, but we only detect the southern half of 
the /-band source in Ha. We explored the reduced data cube 
of the source, and we find (i) there are two clear velocity 
gradients in the cube, reinforcing the presence of two differ- 
ent components and (ii) there is an OH sky-line in the chan- 
nels in which we expect the northern part of the companion, 
explaining the non detection in the maps. In addition, the 
higher velocity dispersion of the MASSIV target occurs in 
the expected overlapping region between both components. 
We conclude that this is a close pair system with r p = 4.2h 
kpc, Av ~ 130 km s _1 and Am, = 0.8. Thus, we classify the 
system as a major merger. 

- 910296626 (Fig. [19]). Major merger. The velocity map 
shows two different components with a separation of r p = 
12.1/i _1 kpc and a relative velocity of Av ~ 165 km s _1 . 
The companion is located toward the north-east. We find 
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Fig. 13. The same as Fig. [3] but for the MASSIV source 020283083 (major merger). The outer contour marks the 2.12 ADU (4<x s ky) 
isophote, while brighter isophotes increase in 1.06 ADU (2cr sky ) steps. [A colour version of this plot is available at the electronic 
edition]. 
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Fig. 14. The same as Fig. [3] but for the MASSIV source 020283830 (minor merger). The outer contour marks the 1.14 ADU 
(3.5cr sky ) isophote, while brighter isophotes increase in 0.75 ADU (2.3cr sky ) steps. [A colour version of this plot is available at the 
electronic edition]. 
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Fig. 15. The same as Fig. [3] but for the MASSIV source 020465775 (major merger). The outer contour marks the 2.92 ADU (9cr sky ) 
isophote, while brighter isophotes increase in 0.97 ADU (3cr sky ) steps. [A colour version of this plot is available at the electronic 
edition]. 



Ami - -0.1 and Am^, = -0.2, so the system is a major 
merger. 

- 910337228 (Fig. [20}. Major merger. The velocity map 
shows two different components with a separation of r p = 
9.5h~ l kpc and a relative velocity of Av ~ 220 km s~'. The 
companion is toward the west and has Am,- = 1 .4. Thus, we 
classify the system as a major merger. 



In summary for this redshift range, we classify 7 of the 1 1 
close pairs candidates as major mergers, and the four other ones 



as minor mergers. This translates to a gas-rich major merger 
fraction of /mm = 0.201!°;°*° at T ia = 1.32. 



4.3. Close pair candidates at 1.5 < z < 1.8 

The weighted mean redshift of the third redshift bin is T r 3 = 
1.54. This is a redshift range where there is no measurement 
of the merger fraction from spectroscopic close pairs yet. We 
identify 3 close pair candidates over 12 galaxies: 
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Fig. 16. The same as Fig. [3] but for the MASSIV source 220376206 (minor merger). The outer contour marks the 3.39 ADU (3<x S k y ) 
isophote, while brighter isophotes increase in 2.26 ADU (2cr sky ) steps. [A colour version of this plot is available at the electronic 
edition]. 
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Fig. 17. The same as Fig. [3] but for the MASSIV source 220544103 (major merger). The outer contour marks the 3.32 ADU (2cr S |j y ) 
isophote, while brighter isophotes increase in 3.32 ADU (2cr s k y ) steps. [A colour version of this plot is available at the electronic 
edition]. 



VVDS-Ultradeep: 910154631 




Velocity field 



Velocity dispersion map 




2 = 1.335 : : 

%nfoni = 4800 s 

1.0 (1.5 0.(1 -0.5 -1.0 -1.5 

A a (") 




A a (") 



Fig. 18. The same as Fig. [3] but for the MASSIV source 910154631 (major merger). The outer contour marks the 1.9 ADU (6cr S |j y ) 
isophote, while brighter isophotes increase in 0.8 ADU (2.5cr s ]q,) steps. [A colour version of this plot is available at the electronic 
edition]. 



- 020116027 (Fig. 121). Major merger. The velocity map 
shows two different components with a separation of r p = 
26.8/z -1 kpc and a relative velocity of Av ~ 100 km s _1 . The 
companion is toward the north-west. The luminosity differ- 
ence is Am,- = 0.7, suggesting a major merger. The difference 
in the K s band is A/Kg, = 0.5, confirming the previous major 
merger classification. 

- 910186191 (Figf22|. Major merger. The velocity map 
shows two different components with a separation of r p = 
12.1 hT x kpc and a relative velocity of Av ~ 450 km s~' . The 
MASSIV target is only detected in 6 pixels because there is 



an OH sky-line in the position of Ha at its redshift. The com- 
panion, located toward the west, is well detected. We find 
Am, = -0.2, this is, the companion is slightly brighter than 
the principal galaxy. However, the MASSIV target is barely 
detected in the K s band, with Am^ - -2.4, suggesting a low 
mass system. This is the only system in which the classifica- 
tion in the two bands is different. Fortunately, both sources 
are VVDS targets, and we have an estimation of their stel- 
lar masses from SED fitting. The difference in stellar mass is 
p ~ 1/3, so we classify the system as a major merger. 
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Fig. 19. The same as Fig. [8] but for the MASSIV source 910296626 (major merger). The outer contour marks the 3.32 ADU (lOcr^) 
isophote, while brighter isophotes increase in 2.56 ADU (So^ky) steps. [A colour version of this plot is available at the electronic 
edition]. 
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Fig. 20. The same as Fig. [8] but for the MASSIV source 910337228 (major merger). The outer contour marks the 1.73 ADU (4-o'^y) 
isophote, while brighter isophotes increase in 1.30 ADU {3cr^y steps. [A colour version of this plot is available at the electronic 
edition]. 
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Fig. 21. The same as Fig. [8] but for the MASSIV source 020116027 (major merger). The principal galaxy is that in the south-east. 
The outer contour marks the 2.02 ADU (3.5<x s ]q,) isophote, while brighter isophotes increase in 1.16 ADU ^cr^y) steps. [A colour 
version of this plot is available at the electronic edition]. 



- 910274060 (Fig. [23j. Major merger. The velocity map is 
consistent with one single component. However, the position 
angle (PA) from the /-band photometry, PA = 105° (North 
has PA = 0° and East has PA = 90°), is nearly perpendicu- 
lar to that from the kinematical modelling, PA = 33°. This 
suggests a complex system, so we performed the GALFIT 
modelling with two sources. We recover well two sources, 
one in the north and the other in the south. The photomet- 
ric PAs of these two sources, provided by the GALFIT fit- 
ting, are now in better agreement with the kinematical one 
(PAi = 75°, PA 2 = 41°), supporting that this is a close pair 



system. The separation between the components is 3. Ah 
kpc, with Av ~ 10 km s _1 , and the luminosity difference is 
Am, = 0.2. Thus, we classify the system as a major merger. 

In summary, we identify the 3 candidates as major mergers. 
This translates to a gas-rich major merger fraction of /mm = 
0.323^q'jq^. Note that in this range our merger candidates have 
r max _ 30/j-i t0 improve the statistics. Applying Eq. O we 
estimate / MM = 0.220!°;^ for r™ x = 20ft- 1 kpc at ^ >3 = 1.54. 

Our results alone, summarized in Tabled suggest a constant 
major merger fraction of /mm ~ 0.21 at 0.9 < z < 1 .8 for r™ ax = 
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Fig. 22. The same as Fig. [3] but for the MASSIV source 910186191 (major merger). The outer contour marks the 2.83 ADU (8<x sky ) 
isophote, while brighter isophotes increase in 2.83 ADU (8a" s k y ) steps. [A colour version of this plot is available at the electronic 
edition]. 
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Fig. 23. The same as Fig. [3] but for the MASSIV source 910274060 (major merger). The outer contour marks the 1.20 ADU (3cr S |j y ) 
isophote, while brighter isophotes increase in 1.00 ADU (2.5o" s k y ) steps. The grey dot in the central figure marks the centre of the 
system given by the kinematical model, while the red/blue dot marks the photometric centre of the principal/companion galaxy 
given by the two-component model of the source. The long bars mark the position angle of the previous models. The short bar 
marks the photometric position angle given by the one-component model of the source. [A colour version of this plot is available at 
the electronic edition]. 



Table 2. Gas-rich major merger fraction and rate of star-forming galaxies at 0.9 < z < 1.8 in the MASSIV sample. 



Zr 


N 




Zr 


log (MJM a ) 


r max 
P 

(h- [ kpc) 


7"mm 
(Gyr) 


/mm 


(Gyr 1 ) 


0.94 < z < 1.06 
1.2 <z< 1.5 
1.5 <z< 1.8 


18 
30 
12 


3 
7 
3 


1.03 

1.32 
1.54 


10.17 
10.57 
10.09 


20 
20 
30 


1.80 
1.37 
2.54 


0.208+°.'! 

0.20 r™ 

U,J "-0.I07 


11 6+° 084 

a 1 47+0.058 
"•i^' -0.037 
n 1 97+O.O79 
-0.042 



20/i 1 kpc close pairs (Fig. l24l>. This merger fraction at z > 1 
is higher by an order of m agnitude than in the local universe 
where f MM ~ 0.01 - 0.03 (IPatton et alJl2007j |Pe Propris et al 



2007 
2010; 



Patton & Atfield 2008} iDomingue etaT]|2009t iDarg et al 



Xu et al.l l2012). This is the first main result of the present 



paper. We compare our major merger fractions with others in the 
literature in Sect. 16. II 



5. The gas-rich major merger rate in MASSIV 

In this section we estimate the gas-rich major merger rate (Rmm), 
defined as the number of mergers per galaxy and Gyr, of star- 
forming galaxies at 0.9 < z < 1.8. We remind here the steps 



to transform a mer ger fraction to a merger rate. Following 



Ide Ravefet 



n a merg« 
aDd2009), 



we define the major merger rate as 



/MM Cm 7"mM' 



(6) 



where the factor C m is the fraction of the observed close 
pairs that finally merge in a typical time scale 7mm- The 
typical merger time scale can be estimated by cosmological 
and /V-body simulations. In our case, we compute the ma- 
jor merger time sca l e from the cosmological simulations of 
Kitzbichler & White (2008), based on the Millennium simula- 
tion (Spri ngel et al.ll2005l) . This major merger time scale refers 
to major mergers (fi > 1/4 in stellar mass), and depends mainly 
on r™ ax and on the stellar mass of the principal galaxy, with 
a weak dependence on redshift in our range of interest (see 
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Ide Ravel et alJ I2009L for details). We measured the median- 
weighted stellar mass from MASSIV sources in each of the 
three redshift bins under study, and estimated the merger time 
scale for these stella r masses. These time scales already include 
the factor Cm (see iPatton & Atfieldl l2008t iBundv etafl l2009t 
iLin et alj|2010t lL6p ez-Saniu an et alj|201 1 ), so we take C m = 1 
in the following. In addi tion. lLope z-Saniu an et al.1 (1201 ll) show 
that the time scales from lKitzbichler & Whi te (2008) are equiv- 
alent to th a t from the A^-body/hydrodynamical simulations by 
lLotz et all d2008l) . However, we stress that these merger time 
scales have an a dditional factor of two uncertainty in their nor- 
malization (e.g., Hopkins et al. 2 lot lLotz et al.ll201 lb . We sum- 
marize the stellar masses, the merger time scales and the gas- 
rich major merger rates in Table [2] As for the merger faction, 
MASSIV data suggests a nearly constant major merger rate at 
0.9 < z < 1.8, R M m ~ 0.12 Gyr 1 (Fig.|25). We study in details 
the evolution of the major merger rate at z < 1 .5 in Sect. 16.21 



6. The redshift evolution of the gas-rich major 
merger fraction and rate up to z ~ 1.5 

In this section we use the MASSIV results at z > 1 to expand the 
study of the gas-rich major merger fraction (Sect. 16.11 ) and rate 
(Sect. 16.21 ) from spectroscopic close pairs to the redshift desert. 
Then, we explore the importance of gas-rich major mergers in 
the assembly of the red sequence since z ~ 1.5 in Sects. l6~3l and 



6.1. The redshift evolution of the gas-rich major merger 
fraction 

In this section we compare the merger fraction from MASSIV 
with those from previous works. Because the merg er fraction 
evolution depe nds on mass (i .e., de Ravel et al. 2009. [201 ll). lu- 
minosity (i.e.. Ide Ravel etal.ll2009l: lLop ez-Saniuan et alj|2010ah 
and colour (i.e.. ILin et al.ll2008b IChou et al.ll2.01 ll) . we focus on 
samples with M* ~ 10 10 ~ la5 M (ISalpeterfl 19551 IMF) to mini- 
mize systematics. In addition, this mass regime is greatly dom- 
inated by gas-rich (wet) merge rs, as those that we observe in 
MASSIV, at least at z > 0.2 (ILin et al.1 120081; Ide Ravel et al l 
l2009HChouetai1l2011h . 

We define the major (/i > 1 /4) merger fraction normalized 
to r™ x = 20/T 1 kpc as 

/20/T 1 kpc\ 095 

/ MM (20, 1/4) = C p FQi) ^ J fmA{r^,n), (7) 

where the factor C p = r p nax /(r p imx - r™") accounts for the miss- 
ing cl ose companions at small radii in those studies with r™ ln > 
(e.g jBell et al.ll2006l) . and the factor F(jj) translates the merger 
fraction for a given ft to the major merger fraction using the /mm 
depe ndence on u de rived bv lLopez-Saniuan et al.l (l201ll 12012) 
and lXuet all d2012l) . /mm 

oc /i 9±0 4 . The search radius depen- 
dence of t he major merger frac t ion, / mm oc '"p 95 , is the empirical 
one from lLopez-Saniuan et alJ d201 ll) . With Eq.UJwe avoid sys- 
tematic differences due to the close pair definition when compar- 
ing differentworks. 

Ide Ravel et al.l (120091) study the major merger fraction of 
M* > 10 9 - 75 M Q (M* ~ 10 10 - 25 Mp) galaxies in VVDS-Deep 
by spectroscopic close pairs, while lLopez-Saniuan et all d201 lh 
provide the major merger fraction of blue (star-forming) galax- 
ies with ~ 10 10 - 55 M in the same sample. In both studies 



0.3 - 



0.05 



0.01 



. 1 1 
/mm oc (1 - 


I I 








• MASSIV 

AT VVDS-Deep 
► DEEP2 
O COSMOS 

t 2MASS/SDSS 



0.0 



0.5 



1.0 



1.5 



Fig. 24. Gas-rich major merger fraction of M* ~ io 10 ~ 10 - 5 M Q 
galaxies as a function of red shift. Circles ar e from this MASSIV 
data set, triangles are fro m de Ravel et al. (2009) and inverted 
triangles are from lLopez-Saniuan et al. ( 201 ll). both in VVDS- 
Deep, right-pointing triangles are from ILin et al. | (120081) in 
DEEP2 redshift survey, pentagons are from lXu et alj d2012h in 
the COSMOS field, and the hexagon is from lXu etlil] 020121) in 
2MASS/SDSS. The solid line is the least-squares fit of a power- 
law function, /mm = 0.0066x(l +z) 3 91 , to the data. The grey area 
marks the 3cr confidence interval in the fit. [A colour version of 
this plot is available at the electronic edition]. 



r max _ 1 oo/j 1 kpc. lLin et aD(l2008l) report the number of com- 
panions of -21 < M B + 1.3z < -19 galaxies (M* ~ 10 10 25 M Q ) 
with I0 h~ l kpc < r p < 30/T 1 kpc in three DEEP2 redshift 
survey dDavis et alj|2003l) fields. Their principal and compan- 
ion sample are the same, so they miss major companions near 
to the selection boundary. Thus, we apply an extra f actor 1.74 
to Eq . [71 to account f or these missing companions (Linetal. 
l20M . IXiT^t^(l20T2l) measure the fraction of galaxies in close 
pairs wi th u > 1 /2.5 in the CO SMOS^ (Cosmological Evolution 
Survey, IScoville et all 120071) and SDSS0 (Sloan Digital Sky 
Survey. IXbazaiian et all 120091) surveys for M* ~ 10 10 2 M Q 
galaxies. We applied a factor 0.5 to pass from their number 
of galaxies in close pairs to the number of close pair systems 
in the sample (C. K. Xu, private communication), and a factor 
F( 1/2.5) = 1.5 + 0.3 to obtain the major merger fraction. All 
these published (gas-rich) major merger fractions are shown as a 
function of redshift in Fig. (l24b . together with the values derived 
from MASSIV. 

We parametrize the redshift evolution of the (gas-rich) major 
merger fraction with a power-law, 



/mm - /mm,o (1 + zf 



(8) 



The least-squares fit to all the data in Fig. [24] yields /mm,o = 
0.0066 + 0.0006 and m = 3.91 + 0.16. We find good agree- 
ment between all works, with the MASSIV point at z ~ 1 being 
higher than expected fr om the fit, but co nsistent within errors 
with the measurement of ILin et al ] d2008l) at that redshift. In the 
next section we show that this difference disappears when the 
stellar mass of the samples is taken into account, emphasizing 
the importance of comparing results from similar parent sam- 
ples. 



4 http://cosmos.astro.caltech.edu 

5 http://www.sdss.org/ 
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Fig. 25. Gas-rich major merger rate of M* ~ 10 M Q galax- 
ies as a function of re dshift. Circles a r e from this MASSIV data 
set, triangles are from de Ravel et al. 1(120091) and inverted trian- 
gles are from lLopez-Saniuan et al. ( 201 lb. both in V VDS-Deep, 
right-pointing triangles are from Lin et al. (2008) in DEEP2 
redshift survey, squares are from [Lopez-Saniuan et al.l (|2009b) 
in GOODS-S f rom morphological criteria, pentagons are from 
XuetalJ d2012l) in the COSMOS field, and the hexagon is from 
Xu et al.1 (120121) in 2MASS/SDSS. The solid line is the least- 
squares fit of a power-law function, Rmm = 0.0048 x (1 + z) 3 ' 95 , 
to the data. The grey area marks the 3cr confidence interval in 
the fit. [A colour version of this plot is available at the electronic 
edition] . 

6.2. The redshift evolution of the gas-rich major merger rate 

We use Eq. ((6]) to translate the original (i.e., without any nor- 
malization in r™ x ) major merger fractions reported in previous 
section into merger rates. We show them in Fig. [25] The good 
agreement between different works is remarkable, reinforcing 
the idea that the merger time scales used account properly for 
the dependence of the merger fraction both on r™ ax and on stel- 
lar mass. 

We als o show the major merger rate from morphological cri- 
teria from [Lopez-Saniuan et alj 32009b). They measure the gas- 
rich merger fraction of ~ 10 10 5 M galaxie s from asymme- 
tries (A) in the GOODS0 dGiavalisco et al.ll2004l) South field and 
take into account the effect of observational errors in z and A, that 
lead to overesti mations in the major merge r fraction by a fac- 
tor of two-three dLopez-Saniuan et alj2009a||b[). using maximum 
likelihood techniques developed in Lopez-Saniuan et al. (2008). 
Other studies find good agreem ent between the asymmetry- 
based major merger rates from Lopez-Sanjuan et al. (2009b) 
and those from close pair statistics (|L6pez-Saniuan et aLll2010at 
lLotz et al.l201 Ulde Ravel et al.l201 lb . confirming the robustness 
of their methodology. 

We also parametrize the redshift evolution of the (gas-rich) 
major merger rate with a power-law, 



R 



MM 



«MM,0(1+Z)". 



(9) 



The least-squares fit to all the data in Fig. [25] yields Rmm.o - 
0.0048 + 0.0003 Gyr 1 and n = 3.95 ± 0.12. The agreement 
between different works points out the importance of comparing 
results from sim ilar parent samples to avoid systematics (see also 
lLotz et al1l201 ll for an extensive discussion on this topic). We 
also point out that the merger rate and the merger fraction of 
star-forming galaxies show a similar evolution with redshift. 



The second main result in this paper is that the major 
merger rate is well described by a power-law function up to 
z ~ 1.5. However, we note that our MASSIV data seem 
to indicate a flattening of the merger rate evolution beyond 
z ~ 1. Previous studies from mor phological criteria (e.g. 



IConselice 20031: Conselice et al.l 120081) and from photometric 
pairs (Rya n et al . 2008) suggest that the power-law parametriza- 
tion is not longer valid at z > 1.5, where a lower merger fraction 
than expected from the low-z evolution is measured, possibly 
indicating a maximum in t he major merger rate at z ~ 2 (e.s 
IConselicell2006t iRvan et"aill2008t iLopez-Saniuanet al. 20091 
Our new measurements agree with this picture and measure- 
ments from spectroscopic close pairs beyond z ~ 1 .5 are needed 
to test the early evolution of the merger fraction. 

The power-law index n = 3.95 + 0.12 is higher than several 
previous measu rements in the literature (e.g., Bridge et al. 2010; 
Lotz et al l20Tlb . as well as our major merger fraction evolution, 
m - 3. 9 1+0. 16. However, our results refer to M* ~ 10 10 " ia5 M o 
star-forming galaxies, and it is known that the merger fraction 
and rate evolve faster for blue, st ar-forming galaxies than for the 
red and global populations (e.g.. [L"in et alj|2008t Ide Ravel et ail 



I20M IChou et alJlfoTl lLopez-Saniuan et alJl2011h" 



6.3. Number of gas-rich mergers since z = 1.5 

We can obtain the average number of gas-rich major mergers per 
star-forming galaxy between Z2 and z\ < Z2 as 



Nmm(Zi,Z2) 



-r 



Rmm dz 
(l+z)H E( Z y 



(10) 



http://www.stsci.edu/science/goods/ 



where E(z) = y&A + ^m0 + z) 3 in a flat universe. Using the 
merger rate parametrization in Eq. d5), we obtain /Vmm(0, 1.5) = 
0.35 + 0.04. Interestingly, half of this merging activity happens 
atz > 1, with/V M M(l, 1-5) = 0.18 +0.02 and /V MM (0, 1) = 0.17 + 
0.02. Because the cosmic time lapse in these redshift intervals is 
1.55 Gyr and 7.7 Gyr, respectively, the average merger activity 
was higher at 1 < z < 1.5 than at 1 < z by a factor of five. 
In the next section we further explore the consequences of this 
very different major merger activity above and below z ~ 1 for 
the assembly of the red sequence. 



6.4. Testing the major merger origin of massive early-type 
galaxies 

The number density of massive (M* > 10 n M Q ) early- 
type galaxies (E/S0, ETGs in the followi ng) have in- 
creas e d with cosmic time since z ~ 3 (e.g., iPozzetti et al.l 
l2010t IBuitrago et all l201ll) . with ETGs being the domi- 
na nt population among massive galaxies only since z ~ 
1 dVergani et al.1 120081: IBuitrago et al.1 1201 it Ivan der Wei et al.l 
120 1 lb Ivan Dokkum et al.1 1201 ll) . Gas-rich major mergers 
have been proposed as an efficient mechanism to trans- 
form star-forming late-type galaxies into red ETGs (e.g., 
Naab et al.l 120061: iRothberg & Josephl l2006allbt iHopkins et al.l 
2008; Romberg & Fischer! 120 1 0t iBournaud et al.l 1201 ll) , so the 
comparison between the observed number density evolution of 
ETGs (petg) and the major merger history of star-forming galax- 
ies imposes important constraints on the role of mergers in 

galaxy evolution. 

In their work, iRobaina et all d2010l) and iMan et al.1 (120121) 
integrate the observed major merger rate over cosmic time and 
predict the evolution of the number density of massive galaxies 
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assuming that all the merger remnants are new massive galax- 
ies. They suggest that major mergers are common enough to ex- 
plain the number density evolution of mass ive galaxies (ETGs 
+ spir als) at z < 1 and z < 3, respectively. lEliche-Moral et alJ 
(2010) model the evolution of the luminosity function backwards 
in time for bright galaxies, selected according to their colours 
(red/blue/total) and their morphologies. They find that the ob- 
served luminosity function evolution can be naturally explained 
by the observed gas-rich and dry major merger rates, and that 
50-60% of massive ETGs in the local universe were formed by 
major mergers at 0.8 < z < 1, with a small number evolution 
since z = 0.8. 

In this section we implement a model to explore the role 
of mergers in the number density evolution of massive ETGs 
since z ~ 1.3. As reference values, we use the number den- 
sities of massive (Af* > 10 1125 Af Q ) ETGs provided by 
Buitrago et al.1 (1201 ll) . They perform a consistent morphological 
study by visual inspection between z — and z = 3, combining 
the SPSS, POWIR (P alomar Observatory Wide-Field Infrared, 
Conselice et aT]j2007h and GNS[] (GOODS NICMOS Survey, 
Conselice et alJuOl lal) surveys. We also use the number den- 
sities of massive spheroidal galaxies with Af* > 10 11 25 M n , se- 
lected by automatic indices in the z COSMOsQ ( L illy et alj2 007 ) 
survey, provided by Poz zetti et al.l d2010h . We show these num- 
ber densities in Fig. [26] 

In our model we assume that, after the final coalescence 
of the merging galaxies, a cosmic time Af = 0.5 Gyr is nec- 
essary for the merger rem nant to be classified as an ETG (see 
lEliche-Moral et aO2010l for a detailed summary of this topic). 
This implies that the new ETGs which appeared between z max 
and z m j n came from the merger activity in the range z e {Z\,Zz), 
where z\ = z m in + Az, zi = Zmax + Az, and Az is the redshift in- 
terval that spans At in each case. Therefore, we take z m ax = 1.3 
as the upper redshift in our model because that implies zi ~ 1.5, 
the redshift limit of the present merger rate study (see Sect. 16.21 ). 

The number density of new ETGs with stellar mass Af* > 
Af*,ii m from gas-rich major merger events in the range z m j n < z < 

2max IS 



Af*, hrn ) = OflMM^MM/LTcdM^dz, (11) 

\J z\ *J -co 

where the different elements in the integral are (Fig. [27j: 

- the global mass function, parametrized with a Schechter 
function 



<&(z,M*) 



4>*(z) I M+ 

m;(z) \m:(z) 



<)(;) 



exp 



M* 

m;(z) 



(12) 



We assume the Schechte r function parameters from 
iPerez-Gonzalez et alJ (120081) and use their pa rametrization 
with redshift provided bv lLopez-Saniuan et al.l (1201 Oal) . 



Iog(0*(z)/Mpc- 3 ) = -2.72 - 0.56(z - 0.5), 



log(M;(z)/M ) = 11.23 



a(z) 



-1.22 



0.13(z-0.5), 
- 0.04(z - 0.5). 



(13) 
(14) 
(15) 



The gas-rich major merger rate, Rmm(z), as measured in 
Sect. 16.21 We assume that it is independent of the stellar 
mass. 



7 http://www. nottingham. ac . uk/astronomy/gns/ 

8 http://www.astro.phys.ethz.ch/zCOSMOS/ 



The merger efficiency function, E(z, M+, M+ t \i m , [i max , /u m i n ). 
This function takes into account the dependence of the 
merger fraction on the mass ratio fi, /mm k H s 
(Lonez-Saniua n et aDl201 ll) . between /u m i a and /i max , as well 
as the extra stellar mass in the final remnant due to the star 
formation during the merger (/ s f). Formally, 



1 



E = 



A /disc i 

(1 - /disc) X (jlln 

i 



■mm 



s 



if M* >M*,i, 
ifM*, 2 <M^ <M i ,, u 
if Mi, < M*, 2 , 



where 



M 



*,lim 



(1 +jUmin)(l +/sf)' 

(1 +/WX1 +/sf)' 
_ M^,i im _ 

M "M*(1+/ Sf ) ' 



(16) 

(17) 
(18) 
(19) 



and /ji sc is the fraction of disc-dominated remnants after 
the merger event (i.e., the fraction of rebuilding discs). We 
assume /disc = in the following (see Sect. 16.4.21 for 
further details). We define the major merger efficiency as 
£mm = E(z,M*,M^ hm ,id m;a = 1,//™,, = 1/4), while the 
minor merger efficiency as E mm = E(z, M+, M* i; m ,/i max = 
l/4,/U m i n = 1/10). As is Sect. 16.11 we assume s = -0.9. To 
estimate /f we use the parametrization of the gas fraction 
as a function of stellar mass and cosmic time provided by 
Rodrigues et al.l d2012l) . 



/gasfe Af*) 



1 



+ Af , 



gas 



1 +(Af*/10 A «) B »' 



(20) 



where t is the cosmic time between redshift z and the present 
and 

A(t) = 9.15 + 0.13?, (21) 
B(t) = 0.5 + 13.36 x exp(-38.02/f). (22) 

In add ition, we assume the prescriptions in Hopki ns"et al.l 
(2009) to estimate the amount of the initial gas mass that 
is transformed into stars during the merger, /burst = /gas(l - 
f m )2n/(l+fi). Finally, 



/f = M burst /Af* 



/gas; 



(23) 



where the factor 2/3 is the efficiency for fi - 1/2 
mergers, which is the typical m ajor merger mass ratio 
dLopez-Saniuan et alj|201lll2012l) . 
- The fraction of late-type galaxies (spirals and irregulars), 
/ltg(z> Af*), this is, the number of potential gas-rich merger 
progenitors over the global population. We parametrize this 
fraction as 

/Lxcfc Mi,) = 0.37 + 0.16z - O.5[log(A/*/A/ ) - 1 1.45] .(24) 



Buitrago et al 



Bernardi et al 



We as sume the redshift dependence find by 
(2011), while the mass dependence from 
(2010) in SDSS. This mass dependence is p resent for 
red/blue galaxies up to z ~ 1.2 (lllbert et al.1 12010). 

Finally, we define the fraction of new massive ETGs due to 
jas-rich (wet) mergers as 



/wet(2mim Zmax) — 



Pwetfemin; Zmaxi 10 Af Q ) 
PETG(Zmin) -pETG^max) 



(25) 
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look 
3 



- back time [Gyr] 
5 6 7 8 

T 




0.0 0.2 0.4 0.6 0. 

z 



1.4 



Fig. 26. Number density evolution of massive (M+ > 



10 



11.25 



M ) ETGs ( E/SO) as a functi on of redshi ft from 



iBuitrago et ail d201ll pentagons) and iPozzetti et all (1201 (X 
squares). The dashed line is the expected number density evo- 
lution due to gas-rich (wet) major mergers from our model. The 
solid line is the expected number density evolution due to wet 
major and dry mergers (both major and minor) from our model. 
Mergers are common enough to drive the number density evo- 
lution of massive ETGs since z ~ 1.3. [A colour version of this 
plot is available at the electronic edition]. 




10.8 11.0 11.2 11.4 11.6 11.8 
log (MJM Q ) 

Fig. 27. Number density distribution from the different func- 
tions in Eq. (fTTT i at z = 1. The solid line is the stellar mass 
function (<P) multiplied by the major merger rate (Rmm)- The 
dashed line includes the major merger efficiency function (Emm)- 
The dotted line includes the fraction of late-type galaxies (/ltg) 
and provides the final function that we integrate over stellar 
mass and cosmic time to obtain p wet . The vertical dashed line 
marks the limiting stellar mass that define massive galaxies, 
M^iim = 10 1L25 M , while the dotted lines mark the stellar 
masses M*,i = 10 1015 M and M* ;2 = 10 1105 M when / sf = 
(see text for details). [A colour version of this plot is available at 
the electronic edition]. 



Our merger model finds /„ e t(0, 1.3) ~ 55%, while this frac- 
tion increases up to / we t(z, 1.3) > 90% when we focus on the 
high-z regime, z > 0.8 (dashed line in Fig. l26l >. This indicates 
that gas-rich major mergers are common enough at z > 0.8 
to fully explain the observed increase in the number density 



of massive ETGs. However, at z < 0.8 only one third of the 
evolution can be accounted by these mergers, / wet (0, 0.8) ~ 
40%. That supports the idea that gas-rich major merging is 
the main proces s involved in the a ssembly of the red sequence 
at z > 1 (e.g.. Illbertet alj|2010t lLopez-Saniuan et all 120105 
lEliche-Moral et alJl2010h . 

In addition to wet mergers, dry mergers can also increase 
the number density of massive ETGs. In this case they promote 
ETGs with M+ < 10 1L25 M to the massive regime. We estimate 
the contribution of both major and minor dry mergers as 



PdrvCZmin, Zmax, ^*,lim) 



-n 



O £ M m/etg dM*dz 



— oo 
<W*.l 



0/C 3 ^mm/ETGdM*dz, (26) 



where R^^ (^ m n?) is th e major (minor) merger rate of mas- 
sive ETGs galaxies from lLopez-Saniuan et all d2012l) . /etg = 
1 - /ltg and we assume ff — in both major and minor merger 
efficiency functions. Note that the integratio n in mass space only 
reach M+ ii m and that the merger rates from lLopez-Saniuan et all 
(2012) include also mixed mergers (ETGs - LTGs). However, 
dry mergers between two already massive ETGs decrease the 
number density petg- We take this into account with the follow- 
ing function, 

Piv(Zmin, Zmax, M +Jim ) = I J $ fi™ 6mm/etG dM*dz 
•Jz\ JM tSm 

Oi^T m G emm/ETcdM+dZi (2?) 

where in this case the merger efficiency function has the form 

[ 0.65, ifM*>M* 3 , 
e = \ 0.65 x (p/n mia ) s if M M < M-t, < M*, 3 , (28) 
0, if M+ < A , 



where 



Mi,,3=fx^Mi, tiim , 
Mi, A = M*j im , 
_ Mi, Xlm 



V = 



M^ 



(29) 
(30) 

(31) 



and the factor 0.65 is the fraction of companions of 
massive galaxies that are already early-type/r ed galaxies 
(lLopez-Saniuan et all 12012b iNewman et all 120121) . As previ- 
ously, we define the major merger efficiency as e M M = 
e(z,M*,M*,i im ,^ max = l,ju min = 1/4), while the minor merger 
efficiency as e mm = e(z,M*,M* > i im ,jt/ max = l/4,fi min = 1/10). 
Analogous to the wet merger case, we define 



/dry (Zmin» Zmax) — 



Pdiy(Zmin,Zmax, W l Zi M ) 



PETG(Zmin) ~ PETG(Zmax) 
p; y (Zmm,Zmax, 10 1L25 M ) - p dry (z mm , Zmax, 10 11 " 25 M ) 

PETG(Zmin) -PETG(Zmax) 



(32) 



to explore the role of dry mergers in massive ETGs assembly 
since z = 1.3. We find / dly (0, 1.3) ~ 35%. As shown in Fig. [26] 
dry mergers are more important at recent cosmic times due to 
the increase in the number density of ETGs, in contrast with the 
diminishing importance of gas-rich mergers. For example, we 
have /dry(0.8, 1 .3) ~ 15%, while / dry (0, 0.8) ~ 40%. 
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The combined effect of gas-rich and dry mergers, / tot = / we t+ 
/dry, is able to explain most of the evolution in petg since z - 1.3, 
with / to t(0, 1 .3) ~ 90%. Thus, our model suggests mergers as the 
main process in the assembly of massive ETGs since z - 1.3. 
Two thirds of the number density evolution is due to gas-rich 
major mergers, while one third is coming from major and minor 
dry mergers. 

The measurement of the merger rate at z > 1 .5 is needed to 
fully constraint the role of gas-rich major mergers in the early as- 
sembly of the red sequence, as well as the possible contribution 
of cold gas accretion in this mass assembly. 

6.4.1 . Fast and slow rotators in the local universe 



The results from the SAURON project (Ide Zeeuw etal] [2002) 
propose a kinematical class ification of morpholo gical ETGs 
into fast and sl ow rotators dEmsellem et al.l 2007). Recently, 
the ATLAS 3D dCappellari et all I201 lh survey has observed a 
representative sample of 260 nearby ETGs, finding that the 
fraction of slow rotat ors increases with the dynamical mass 
(Ems ellem et alJl201 lh. T he cosmological simulation analysed 
by Khochfaret al. (2011) suggests that the main difference be- 
tween fast and slow rotators is their average number of major 
mergers, with fast rotators having undergone one major merger, 
while slow rotators have undergone two. We use this fact and the 
merging model developed in the previous section, to predict the 
fraction of slow rotators in the local universe, / s i ow . We simply 
assume that wet major mergers produce fast rotators (first major 
merger event), dry major mergers produce slow rotators (second 
major merger event), and dry minor mergers do not change the 
kinematical state of ETGs. In addition, we took all the ETGs at 
Z = 1.3 as fast rotators. 

With the previous assumptions, we expect / s i ow ~ 55%, 
in good agreement with the ATLAS 3D result of / s i ow ~ 47 - 
75% for ETGs with a dynamical mass M dyn > 10 1125 M 
(Ems eilem et alj|201 lh . Thus, our model also reproduces the lo- 
cal fraction of slow/fast rotators of massive ETGs galaxies, rein- 
forcing the conclusions of the previous section. 



6.4.2. Uncertainties in the model 

The model presented in previous sections has set most of its pa- 
rameters from observational results. However, there are two un- 
constrained parameters that could affect our conclusions. The 
first parameter is /dj sc , the fraction of gas-rich major mergers 
whose remnant rebuild a disc component and do not contribute 
to the increase in petg- We as sumed /disc = 0. and in the fol- 
lowing we justify this selection. Hopki ns et alj (T2009) find that a 
high gas fraction prevents the destruction of the disc component 
after a major merger in their /V-b ody/hydrodynamical simula- 
tions (but see Bour naud et al.l20Tll for a different point of view). 
These simulations suggest that disc rebuild could be an efficient 
process when the gas fraction is / gas > 50%. However, we have 
focused in the massive end of the galaxy population, where the 
gas fractions are lower. Thanks to Eq. ( 1201 ). we can estimate the 
gas fraction of our gas-rich mergers. We find / gas < 30%, justify- 
ing the assumed /di sc = 0. At lower masses than explored in the 
present section, the gas fraction is higher, and disc rebuild could 
be an impor tant process in the formation of bulge-dominate d 
spirals (e.g.. iHuertas-Companv et al] 12010b iPuech et"aTll2012h . 
However, we cannot discard positive values of /d; sc for massive 
galaxies, as we see below. 



The second parameter is the assumed merger time scale, 
which typically has a factor of t wo u ncertainty in their 
normalization (e.g., Hopkins et al.l 120101) . The 7mm from 
iKitzbichler & White! d2008l) are typically lo n ger than others 



Ritzbicnier &l wmte (ZUU8) are typically longer man otners 
in the literature (e.g.. iPatton & Atfieldl 120081; Lin et al l 1201 Ol) 
or similar t o those from /V-body/hydrodynamical simulations 
(lLotz et al.l l2010b| |a|). Thus, we expect, if anything, a shorter 
7mm i which implies a larger number density of ETGs due to 
mergers. Nevertheless, the good description of the petg evolu- 
tion with our merger model strongly suggests that mergers are 
indeed the main process in massive ETGs assembly. Thus, a 
lower 7mm (i.e., a higher merger rate that translates to an ex- 
cess of ETGs) could be compensated by a positive value of /di sc , 
that reduces the number density of ETGs due to mergers. 



7. Summary and conclusions 

Using SINFONI/VLT 3D spectroscopy, we have been able to 
measure, for the first time with spectroscopically-confirmed 
close pairs, the gas-rich major merger fraction and merger rate at 
around the peak in star formation activity at 0.9 < z < 1.8, from 
the MASSIV sample of star-forming galaxies with a stellar mass 
range M* = 10 9 - 10 11 M . In this redshift range we identify 20 
close pairs, and classify 13 as major mergers based on a separa- 
tion r p < 20/i~' - 30/i~' kpc and a relative velocity Av < 500 km 



20.8ig; g -% 



We find that the gas-rich major merger fraction is high, 
' 5 g 2 %, 20.1+®°% and 22.0^ 3 3 7 % for r p < 20/z 1 kpc 
close pairs in redshift ranges z = [0.94, 1.06], [1.2, 1.5) and 
[1 .5, 1 .8), respectively. When compared to measurements at red- 
shifts z < 1, the evolution of the (gas-rich) merger fraction can 
be parametrized as /mm = 0.0066 x(l+z)" ! with m = 3.91+0.16. 
We note that the evolution between z — 1 and z ~ 1 .5 seems to 
flatten out compared to lower redshifts. 

The merger rate has been derived using merger time scales 
from the literature: we find that the gas-rich merger rate is 



1 16 +0.084 



Gyr 1 , 0.147+",$* Gyr 1 and 0.127^;^ Gyr 



at 



http://www-astro.physics.ox.ac.uk/atlas3d/ 



z = 1.03, 1.32 and 1.54, respectively, for merger time scales of 
Tum ~ 1 .5 Gyr. We then find that the (gas-rich) merger rate evo- 
lution for galaxies with stellar mass M+ = io 1() ~ 10 - 5 M Q over 
z = [0, 1.5] scales as (1 + z) n withn = 3.95 + 0.12. 

Using these measurements, we develop a simple model to 
estimate the contribution of gas-rich major mergers to the growth 
of galaxies on the red sequence. We infer that ~ 35% of the star- 
forming galaxies with stellar mass M* = 10 10 - 10 10 ' 5 M Q have 
undergone a major merger since z ~ 1.5. The number of major 
merger events was about 5 times higher at 1 < z < 1 .5 compared 
to z < 1 . Assuming that each gas-rich major merger produces 
a new early-type galaxy, we infer that the number of gas-rich 
major mergers is large enough at z > 1 to explain the increase 
in the number density of massive ETGs, supporting a picture 
where gas-rich (wet) merging is the main process building-up 
the red sequence. While gas-rich mergers become rarer towards 
lower redshifts, the number of dry mergers is steadily increasing, 
and the combination of these two processes accounts for most, 
if not all, of the increase in the number density of massive ETGs 
since z ~ 1.3. Two-thirds of this number density evolution is 
due to wet major mergers, while one-third is coming from major 
and minor dry mergers. These results add further evidence to a 
picture where merging is a major process driving mass assembly 
into the red sequence galaxies. 

We note that minor merging is definitely present in the 
MASSIV sample (see Sect. |4j. However, due to incompleteness 
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in detecting these faint companions, we are not able to assess a 
minor merger rate at these epochs from our data. In the global 
picture of red sequence assembly, we emphasize that a simple 
extrapolation of the minor m erger rate measured up to z ~ 1 by 
iLopez-San iuan et al] d201 ll) . would indicate that from z ~ 1.5 to 
the present, minor mergers with 1/10 < /i < 1/4 are not com- 
mon enough to significantly move spiral galaxies into the red 
sequence. 

To get a complete picture of the life-time effect of ma- 
jor merging on massive galaxies, accurate measurements of the 
merger fraction and merger rate are needed beyond z ~ 2. 
Spectroscopic surveys will remain an important element to pro- 
vide secure identification of merging systems at these early 
epochs. 
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N22h, P hot(Zr)/"22h(l, 0) is the fraction of sources with a photomet- 
ric redshift in the redshift range z r over the total population of 
sources with flag = 1 and in the 22h field. 



(A.l) 



Appendix A: Spectroscopic Success Rate in 
VVDS-Wide fields 

We computed the Spectroscopic Success Rate (SSR) in VVDS- 
Wide 22h and 14h f ields for a given re dshift range by following 
the prescriptions in Ilbe rt et alJd2006h : we compare the number 
of galaxies with flag = 4, 3 and 2 over the total number of galax- 
ies (those with flag = 4, 3, 2, 1 and 0). Because flag = 1 galaxies 
have a 50% reliability and no redshift information is available 
for flag = galaxies, we took advantage of the latest photomet- 
ric redshifts from CFHTLS survey to estimate the number of 
galaxies with flag = 1 and that belong to the redhsift range of 
interest. With the previous steps we estimate the SSR in the 22h 
field as 

CCD , s N 2 2h, S pec(Zr) 

SSR 2 2h(Zr) = — — 77 — , 

JV22h,specUr) + ^V22h,photUrJ 

where A^h.specfe) is the number of galaxies with flag = 4, 3 
and 2 in the redshift range z r , and A^h.photfa) is the number of 
galaxies with flag = 1 and 0, and with a photometric redshift in 
the same redshift range. As the available imaging data in the 14h 
eld is not as deep as the CFHTLS, we show in the following that 
we can assume that the SSR in the 14h field follows the same 
distribution as in the 22h field, for which CFHTLS photometry 
is available. 

First we checked the properties of the stars in these two 
VVDS-Wide fields. The fraction of stars in the 22h (14h) field is 
66% (5 1 %) for flag = 4 sources, 38% (26%) for flag = 3 sources, 
24% (23%) for flag = 2 sources, and 11% (12%) for flag = 1 
sources. So the fraction of stars is similar in both fields for flag 
= 2 and 1 sources, while higher in the 22h field for flag = 4 and 
3. We checked that the normalized distributions of the VVDS- 
Wide stars as a function of their observed /AB _ band magnitude 
in both fields are similar for each flag. We note that stars with 
high confidence flag are brighter than the low confidence ones, 
as expected. The fraction of bright stars (/ab < 21) in the 22h 
(14h) field is 75% (67%) for flag = 4 stars, 49% (43%) for flag 
= 3 stars, 30% (29%) for flag = 2 stars, and 19% (16%) for flag 
= 1 stars. These fractions suggest that there is a higher density 
of bright stars in the 22h area than in the 14h one, that translates 
in a higher fraction of stars for flag = 4 and 3 sources, while 
faint stars have similar densities in both fields, leading to a sim- 
ilar fraction of stars with flag = 2 and 1. Because of this, we 
assumed that the fraction of stars among flag = sources, that 
we estimate in the 22h field from the CFHTLS photometry, is 
similar in 22h and 14h fields. 

The distribution of galaxies at z > 0.9, the redshift range 
in which we are interested on, are also similar in both fields. 
Because of this, we assume that the photometric distribution of 
galaxies with flag = 1 and is similar in both fields, and we use 
those in 22h field to estimate A^i4h, p hot and the SSR in the 14h 
field: 



SSRuhiZr) 



iVl4h,spec(Zr) 



Nl4h,spec(Zr) + /22hfe) X n I4h (l, 0) 



(A.2) 



where «i4h(l,0) is the total number of sources (galaxies and 
stars) with flag = 1 and in the 14h field, and /22/ife) 



Appendix B: Completeness of the close pair sample 

Throughout the present paper we assume that the close pair 
systems detected in the MASSIV sample are representative for 
galaxies with /ab ^ 23.9. The MASSIV sample comprises 
galaxies fainter than this luminosity limit, and in this appendix 
we justify the boundary applied in our analysis. 

The VVDS parent samples of the MAS SIV samp le are ran- 
domly selected in the /ab band dLe Fevre et al.l l2005). Because 
of this, we study the distribution of MASSIV sources as a func- 
tion of the /AB-band magnitude to obtain clues about the com- 
pleteness of the MASSIV sample in our close pair study. We 
show the cumulative distribution of all MASSIV sources and 
of those with a close companion (both major and minor) in 
the top panel of Fig. IB. II Obviously, we are able to detect sin- 
gle MASSIV galaxies at fainter magnitudes than the close pairs 
(24.4 vs 23.8) because in the latter case we have to detect both 
the principal source and the companion galaxy, which is usually 
fainter. 

Lets assume for a moment that the detection curve of the 
MASSIV sources is a step function that is 1 for /ab < /AB.iim 
and for fainter magnitudes. In this ideal case, /AB.iim is defined 
by the fainter galaxy (close pair) detected. Thus, we should only 
be able to detect close pairs at /ab < 23.8, and the measured 
merger fraction in the total MASSIV sample (/ab ^ 24.4) will be 
lower than the real merger fraction because of the missing faint 
close pairs with /ab > 23.8. Of course, the detection curve in the 
MASSIV sample is more complicated than a step function, and 
could depend on redshift, geometry, luminosity, etc. Instead of 
trying to estimate the detection curve for MASSIV sources and 
for close companions to recover statistically the missing close 
pairs (see Patton et al. 2 0001 : Ide Ravel et al.ll2009l for examples 
of this kind of correction), we define the /AB~band magnitude up 
to which the detected close pairs are representative of the total 
MASSIV sample, named /AB,comp> an d in our study we only keep 
those sources brighter than /AB.comp to ensure reliable merger 
fractions. 

The MASSIV sample is a representative subsample of the 
global star-forming population dContini et al.ll2012l) . and we ex- 
pect the MASSIV galaxy pairs to be also a random sample of 
this global population. Therefore, the distributions in the /ab 
band of the total and the close pair MASSIV sources should 
be similar when the close pairs detection was slightly affected 
by incompleteness issues. Following this idea, we performed 
a Kolmogorov-Smirnov (KS) test over the total and close pair 
MASSIV galaxies with /ab ^ ^AB.iim, and explored different val- 
ues of /AB.iim, from 23.0 to 24.5. Then, the completeness magni- 
tude /AB.comp was defined by the minimum in the KS estimator 
D, i.e., where both distributions have the lower probability to be 
different. This exercise provides the curve in the bottom panel 
of Fig. IO that states /AB.comp — 23.9. We repeated this pro- 
cedure with the major merger sample, and we obtain a similar 
distribution of D values, reinforcing our choice. 

In summary, in the estimation of the merger fraction we only 
use those MASSIV galaxies with /ab ^ ^AB.comp = 23.9. This 
ensures that the close pair sample is a random subsample of the 
total MASSIV one, and we avoid any bias related with the shal- 
lower detection curve of close pairs compared with that of the 
total MASSIV sample. 
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23.2 23.4 23.6 23.8 24.0 24.2 24.4 
^AB 

Fig. B.l. Top: Cumulative distribution in the apparent 7ab mag- 
nitude of the total MASSIV sample (green line), the close pair 
sample (red line) and the total MASSIV sample with I AB < 23.9 
(blue line). Middle: The same than before, but the red line 
shows the distribution of the major merger sample. Bottom: 
Kolmogorov-Smirnov estimator D of the total MASSIV sample 
distribution against the close pair (blue solid line) and the ma- 
jor merger (red dashed line) distributions. We computed D for 
those galaxies brighter than a given 7ab magnitude. The solid 
(dashed) horizontal line marks the minimum in the D distribu- 
tion for close pairs (major mergers). The vertical solid line marks 
the 7ab magnitude at the minimum value of D, /AB.comp = 23.9. 
[A colour version of this plot is available at the electronic edi- 
tion]. 
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